Emergence of supersymmetry at a critical point of a lattice model 
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Supersymmetry is a symmetry between a boson and a fermion. Although there is no apparent 
supersymmetry in nature, its mathematical consistency and appealing property have led many 
people to believe that supersymmetry may exist in nature in the form of a spontaneously broken 
symmetry. In this paper, we explore an alternative possibility by which supersymmetry is realized in 
nature, that is, supersymmetry dynamically emerges in the low energy limit of a non-supersymmetric 
condensed matter system. We propose a 2+lD lattice model which exhibits an emergent space-time 
supersymmetry at a quantum critical point. It is shown that there is only one relevant perturbation 
at the supersymmetric critical point in the e-expansion and the critical theory is the two copies of 
the Wess-Zumino theory with four supercharges. Exact critical exponents are predicted. 
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I. INTRODUCTION 

Poincarc invariance is the underlying space-time sym- 
metry of relativistic quantum field theories. There are 
only two mathematically consistent ways of extending the 
symmetry in nontrivial ways due to a no-go theorem 
One is a conformal symmetry and the other, a supersym- 
metry. A conformal symmetry combines the Poincare 
invariance with scale invariance. It is realized in the 
long-distance limit of a massless theory at which all fi- 
nite length scales are scaled out. A supersymmetry is a 
symmetry between a boson and a fermion. Generators of 
supersymmetry Qa, which are called supercharges, are 
spinors and they satisfy the commutation relations, 



(1) 



where a, (3 are spinor indices, P^^ are constants, and 
is the energy-momentum operator. Because super- 
charges are fermionic operators, a boson is transformed 
into a fermion (and vice versa) under supersymmetry 
transformations. Therefore the number of bosonic modes 
is equal to the number of fermionic modes in supersym- 
metric theories. The second commutation relation in Eq. 
^ implies [P^P^ , Qa] = and masses of supersymmet- 
ric partners are identical. Since bosons and fermions 
contribute to quantum effective actions with quantum 
corrections of the opposite signs and they have same 
masses, the effects of quantum fluctuations of bosons and 
fermions are canceled with each other in supersymmet- 
ric theories. Because of this, quantum corrections are 
highly constrained by kinematics. If there are enough 
supersymmetries, there is no quantum correction at all 
(non-renormalization) for some quantities. Due to the 
non-renormalization property, supersymmetry has been 
proposed as a way of stabilizing the hierarchy of vastly 
different mass scales present in the standard model. It 
may also play an important role in the unification of the 
gauge interactions. On the other hand, many supersym- 
metric theories have been studied as toy models where 



a supersymmetry enables one to understand strong cou- 
pling physics rigorously [2i, |3|]. 

Despite the unique mathematical consistency and 
beautiful properties of supersymmetry, nature does not 
exhibit supersymmetry at low energy scales. If nature is 
supersymmetric, it should be spontaneously broken. If 
that is the case, supersymmetry will become manifest at 
a high energy scale. An alternative way of finding su- 
persymmetry in nature may be to go to a low energy 
in condensed matter systems, relying on the principle of 
emergence. 

A new symmetry can emerge in the low-energy limit al- 
though a microscopic model does not respect the symme- 
try. For example, a low-energy effective theory can have 
the full Poincare invariance although the underlying lat- 
tice explicitly breaks rotational symmetry in a condensed 
matter system. Even a gauge symmetry[j| and a general 
covariancej^ Q can be emergent. The emergence of a 
new symmetry can be a characteristic of a new state of 
many-body systems 0- Searching for new states of mat- 
ter in condensed matter systems is becoming an impor- 
tant research avenue as new materials which can not be 
understood in conventional theories are synthesized [1, 
and highly controllable correlated many-body systems 
can be fabricated in cold atom systems [Hi • Therefore 
it would be of interest to find a condensed matter system 
which shows an emergent supersymmetry. 

It has been suggested that supersymmetry can 
emerge in the low-energy limit of a non-supersymmetric 
theory [Tsl. [T^. In 1-l-lD, supersymmetry emerges at 
the tricritical point of the dilute Ising model [l5[. Emer- 
gent supersymmetries play important roles in realizing 
lattice versions of supersymmetric field theories in vari- 
ous dimensions flGj. For example, the A/" = 1 3+lD su- 
per Yang-Mills theory can emerge without an underlying 
supersymmetry although the notion of the emergent su- 
persymmetry is rather obscure in this case due to the 
opening of a mass gap caused by confinement. Super- 
symmetric field theories can be also realized in lattices 
by fine tunings of bare parameters [l7j or by a dynam- 
ical mechanism where some supersymmetries which are 
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already present in lattices guarantee the emergence of 
continuum supersymmetries without fine tuning (I^. In 
this paper, we construct a 2+ ID lattice model where su- 
persymmetry may dynamically emerge at a critical point 
without any lattice supersymmetry. 

It is noted that non-relativistic supersymmetries have 
been considered in condensed matter systems [l9l [20| . 
In non-relativistic systems, supercharges are scalars (not 
spinors) and the anti-commutator of supercharges gener- 
ates only energy (not momentum), that is, {Q, Q^} = H, 
where iJ is a Hamiltonian. In such systems, supersym- 
metries play the roles as in 0-l-lD quantum mechanical 
systems [2l|. The present work concerns an emergence of 
a full space-time supersymmetry in a 2-l-lD relativistic 
system where the relativity is also emergent out of a non- 
relativistic microscopic system. In this case, the algebra 
of supercharges generate the translations in both time 
and space through the energy-momentum operator as in 
Eq. (HD. 

II. MICROSCOPIC MODEL AND 
LOW-ENERGY EFFECTIVE THEORY 

The microscopic system is a mixture of fermions and 
bosons. The Hamiltonian is composed of three parts, 

H = Hf + Hb + Hfb, (2) 

where 

= -tf E (///^ + '^•c-), 

<i,3> 

<i,.i> I 
Hfb = ^0 ^ e*"^ (/z+bi //-bi + //-b2 //+b2 + 

+ //-bi+b2//+bi-b2) + h.c. (3) 

Here Hf describes spinless fermions with nearest neigh- 
bor hopping on the honeycomb lattice; Ht describes 
bosons with nearest neighbor hopping and an on-site 
repulsion on the triangular lattice which is dual to the 
honeycomb lattice; and Hfb couples the fermions and 
bosons. The lattice structure is shown in Fig. [T] (a), fi 
is the fermion annihilation operator and e~'^^, the low- 
ering operator of nj which is conjugate to the angular 
variable 9j. i,j and I, J are site indices for the honey- 
comb and triangular lattices, respectively. i/,tb > are 
the hopping energies for the fermions and bosons, respec- 
tively and U is the on-site boson repulsion energy. Note 
that the boson hopping is frustrated. This will play a 
crucial role for the emergent supersymmetry as will be 
shown later, bi and b2 are vectors which connect a site 
on the triangular lattice to the neighboring honeycomb 
lattice sites as is shown in Fig. [T] (a), ho is the pairing 
interaction strength associated with the process where 
two fermions in the f-wave channel around a hexagon are 



paired and become a boson at the center of the hexagon, 
and vice versa. In this sense, the boson can be regarded 
as a Cooper pair made of two spinless fermions in the f- 
wave wavefunction as is shown in Fig. [T] (b) . This model 
has a global U(l) symmetry under which the fields trans- 
form as fi -> f^e^'P and e~''^' -> e~'^'e'^^. 




(b) (c) 



FIG. 1: (a) The lattice structure in the real space. Fermions 
are defined on the honeycomb lattice and the bosons, on the 
dual triangular lattice, ai, a2 are the lattice vectors with 
length a, and bi, b2, two independent vectors which connect 
a site on the triangular lattice to the nearest neighbor sites 
on the honeycomb lattice, (b) The phases of a fermion pair 
in the real space, (c) The first Brillouin zone in the momen- 
tum space. A and B indicate two inequivalent points with 
momenta = ir(|> "7f ) S'ncl ks = ir(|''^) where the low 
energy modes are located, tpi, cj>2 are located at kA and i/'2, 
at kfl. 

First, we identify low-energy modes of the 
fermions and bosons in the absence of the coupling 
Hfb- At zero chemical potential, the fermions are 
half-filled and their energy spectrum is given by 

= ±tf^J(l + cos ki + cos fc2)^ + (sin ki — sin 
where fci — akx, k2 — a +^^v -with a, the lattice 
spacing. There exist two Fermi points at — ^(^, 

and ks = shown in Fig. [1] (c). Since the 

energy dispersion is linear near the Fermi points, the 
low-energy excitations are described by two Dirac 



3 



fermions. 



modes. The resulting interaction Lagrangian is 



(4) 



i=l 



Here ^pi {1P2) denotes the two-component complex 
fermion at momentum k^i (k^)- 9^ = {dr,dx,dy) are 
the derivatives in imaginary time and the spatial direc- 
tions. 7o = (T3, 71 = (Ti and 72 = (T2 with (7^, the Pauli 
matrices. 4'n = ~*'0i7o and c/ '--^ f/a is the fermi veloc- 
ity. 

To obtain a low-energy theory for the bosons, we intro- 
duce a soft boson field $/ = |<I>/|e~*^^ and the potential 
— U2|$P 4- 'ii4|<i>|'* which gives a finite amplitude 
to the soft boson field. In energy-momentum space, the 
boson action becomes 



Sh = 



dk ( ±k 



U2 



+ U4 dkidk2dq $ 



k2 



(5) 



where J dk = J is the energy-momentum in- 

tegration and = 2tb [cos fci + cos/c2 + cos(fci -I- ^2)]. 
Since the boson hopping has the wrong sign, k = (0, 0) 
is not the minimum of e^; rather, two minima occur 
at k^ and k^ where the nodal points of the fermions 
are located. Therefore we have two low-energy boson 
modes. We introduce (/>! and 02 to represent the low- 
energy modes near the kg and k^ points, respectively. 
Note that the (f>i {(f>2) boson carries the same momentum 
as the '4>2 (V'l) fermion. With this convention, we will 
see that only those bosons and fermions which carry the 
same index (n — 1,2) interact with each other at low 
energies. Expanding near the two minima, we obtain 
the effective Lagrangian for the low-energy bosons. 



1=1 



+AiAtEl0„|4 + A2«:|0in</)2|', 



(6) 



where Cf, ^ \/IbU a is the boson velocity which is in gen- 
eral different from the fermion velocity Cf. Although 
both the fermions and bosons have the 'relativistic' en- 
ergy spectrums, there is no Lorentz symmetry if the ve- 
locities are different. The Lorentz symmetry requires the 
velocities of all massless particles to be identical. As will 
be shown later, the Lorentz symmetry will emerge in the 
low energy limit through quantum corrections, m is the 
boson mass and the coupling constants Ai, A2 are made 
dimensionless by introducing a mass scale k. Note that 
momentum conservation does not allow an interaction 
such as (j)24>24'i'Pi- 

We can obtain the interaction between the low- 
energy fermions and bosons by rewriting Hfb in energy- 
momentum space and keeping only the low-energy 



where e is the 2x2 antisymmetric matrix with £12 = 
—£21 = 1. Terms like (t)2''P'[eipi or 02^fe'(/'2 are not al- 
lowed because they do not satisfy momentum conserva- 
tion. 



III. RENORMALIZATION GROUP ANALYSIS 

Now we perform a one-loop renormalization group 
(RG) analysis in 4 — e dimensions for the low-energy ef- 
fective theory given by 



C — Cf + Cb + Cfb- 



(8) 



We use the dimensional regularization scheme where the 
number of fermion components and the traces of gamma 
matrices are fixed[23i|. Maintaining the same number of 
fermionic and bosonic modes in 4 — e dimension is impor- 
tant because supersymmetry requires that the number 
of modes is the same for the bosons and fermions. If 
there is a gauge symmetry, more sophisticated regular- 
ization scheme is necessary to preserve both gauge sym- 
metry and supersymmetry because the number of com- 
ponents of gauge boson should depend on the dimension 
of space-time 2J| . Since there is no gauge symmetry in 
the present model, the simple dimensional regularization 
scheme can maintain supersymmetry [25|. Of course, the 
present model has no supersymmetry. The point is that 
it is convenient to use a regularization scheme which can 
maintain supersymmetry in probing an emergent super- 
symmetry in the low energy limit. 

In the e-expansion, the above Lagrangian contains all 
the relevant and marginal terms. A four fermion inter- 
action has the scaling dimension D — 6 — 2e + O(e^) and 
can be ignored for a small e. In the following, we do not 
consider the four fermion interaction. However, in prin- 
ciple, the four fermion interaction can become important 
in 2 + ID due to a strong interaction, in which case one 
needs to tune a microscopic four fermion interaction term 
to reach the fixed points we will discuss in the following. 
The boson mass is always a relevant perturbation and 
we tune it to zero in order to examine the RG flow of 
the other couplings in the massless subspace. At the one 
loop-level, there are 8 diagrams which are shown in Fig. 
[21 Each diagram contributes to the quantum effective 
action as follows, 



4/i2 



^ i^n folodr + C//l E lidi tpn, 
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(a) 



(b) 




(d) 



(g) 
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(e) 
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'v, 

(h) 
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FIG. 2: One-loop diagrams, (a) and (b) are the self-energy 
corrections of fermion and boson respectively, (c), (d), (e) 
and (f) contribute to the vertex correction of Ai and (f), (g) 
and (h), to the vertex correction of A2. 
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(4.)%gl^"l'' 



I6A1A2 



^^2 I i |2 



(4^)2c2, 



(47r)2c2e 



I0l|1'/'2|^ 



9X2 



(47r)2c^e 



(9) 



where /o = and /i = ^^^^ with a = Cb/c/. 

From the renormahzed quantum effective action, the beta 
functions are obtained to be 



dh 



16c3 



(4^C/)2 



Cb{Cf + CbY 



dXi 
~dl 



1 / 20A? + A2 S/i^Ai 16/i'' 



dA2 

dl 
dcb 

la 



eAo 



(47r)2 
1 



4Ai 



I6A1A2 



8/l2A2 



32h'^Cficb - Cf) 
3(47r)2cb(cf, + c/)2' 

2/.2cb(cg-4) 

(47rc6C/)2 



(10) 



where the scaling parameter I increases in the infrared. 

There are two solutions for /3h = 0. One is the unsta- 
ble solution with h — and the other, the stable one 
with a finite h. At ft, = 0, the bosons and fermions 
are decoupled. The fermion system consists of non- 
interacting Dirac fermions. The RG flow of the boson 
couplings in the subspace of m = /i = is shown in 
Fig. [21 In the subspace of to = ft. = 0, there are three 
fixed points, that is, the Gaussian (GA) fixed point with 
(ft*,At,A^) = (0,0,0), the Wilson-Fisher (WF) fixed 

point with {h*,Xl,Xi.) = (0, ^^"'^f^ O), and the 0(4) 

fixed point with (ft*,A|,A5) = [0, ^^^^). Be- 

cause the linear term of the beta function at the 0(4) 
fixed point accidentally vanishes along a direction, at 
higher order in e there occurs a stable fourth fixed point 
between the WF fixed point and the 0(4) fixed point 
However, the fixed points in the to = ft = plane are all 
unstable because the pairing interaction ft is relevant. 
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FIG. 3: The schematic RG flow of the bosonic couplings in 
the subspace of m = ft = 0. 

Once we turn on the pairing interaction, h flows to 

a finite value with Ii^ — t^'^y) ^ (2 -\ 7— r^— 1 and 

the boson and fermion velocities begin to flow as can be 
seen from the last two equations in Eq. (fTO|) . Because 
the pairing interaction mixes the velocities of the boson 
and fermion, the difference of the velocities exponentially 
flows to zero in the low-energy limit as is shown in Fig. 
[3] (a). The line of Cb = c/ is critical and the value of 
the velocity in the infrared limit is a non-universal value. 
This implies that the bosons and fermions have the same 
energy dispersion and Lorentz symmetry emerges at low 
energies due to quantum fluctuations. Now we consider 
the flow of ft, Ai and A2 with a fixed Cb = Cf = c. The 
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(b) 

FIG. 4: The schematic RG flows of (a) the velocities with 
h^Q and (b) Ai, A2 and h in the subspace of m = 0. In (b), 
the sohd hnes represent the flow in the plane of (/i,Ai) and 
the dashed lines, the flow outside the plane. 



fermions are mixed under the transformation. Since the 
two sets of modes (0„, V'n) are decoupled at the critical 
point, the supersymmetry transformations are indepen- 
dent for n = 1 and 2. That is why the spinor has the 
index n. Here S^^ipn 7^ ~*('^C„V'n)Wo because we are us- 
ing the imaginary time formalism. The supersymmetry 
leads to conserved supercurrents, 

h —T 

= V'„7M#n + *^C^'/'Ie7M- (13) 

Here the spinor indices in and jj^ are suppressed. 
For each n, there are four independent supercharges, 
= Jdx^J^^ and C = / dx^Xa in 2+lD. The su- 
percharges satisfy the commutation relations like Eq. ([T]). 
This corresponds to A/" = 2 supersymmetry in each sec- 
tor of n, that is, there are twice as much supercharges 
as the minimum number of supercharges in 2+lD. The 
resulting super-Poincare invariance is emergent because 
the microscopic model has neither Poincare symmetry 
nor supersymmetry. Each set made of one complex bo- 
son and one two-component fermion forms a chiral mul- 
tiplet of the supersymmetry. This critical theory is the 
Af = 2 Wess-Zumino theory[2^ with two copies of chiral 
multiplets[23|. The critical exponents calculated in the 
one-loop level 

V4> = Vi' = e/3 (14) 



RG flow is displayed in Fig. [3] (b). In the following we 
will use units where c = 1. With a nonzero h, the system 
flows to a stable fixed point, 



12 



12 



,0). 



(11) 



The nonzero h is a consequence of strong pairing fluc- 
tuations at the critical point. This is crucial in obtain- 
ing Lorentz symmetry and supersymmetry as will be dis- 
cussed in the foUowings. At the fixed point, A2 vanishes 
and there is no coupling between the two sets of low- 
energy modes, [ipmi^n) with n = 1,2. Physically, this 
implies that the bose condensates which carry the differ- 
ent momenta k^i and develop independently in the 
condensed phase. At the critical point, Ai = /i^ and the 
theory becomes invariant under the transformation, 



n Sn ^ 7 



(12) 



match those of the Wess-Zumino theory in the leading 
order of e. In the supersymmetric theory, the one-loop 
result is exact as will be explained later. It is of note 
that the critical exponent is independent of regularization 
scheme although the values of the couplings depend on 
regularization scheme. 

The schematic phase diagram of the Hamiltonian ^ 
for a generic value of ho is shown in Fig. [5] in the pa- 
rameter space of tb/U. There is a second order phase 
transition between the normal phase for small tt/U and 
the bose condensed phase for large ti,/U. In the nor- 
mal phase, the fermions are gapless while a gap opens 
in the bose condensed phase. The critical point is de- 
scribed by the supersymmetric Wess-Zumino theory al- 
though both the normal and the bose condensed phases 
are non-supersymmetric. 



Normal phase 



Bose condensed phase 



u/u 



Supersymmetric 
critical point 



where <^ — ^^d^ and ^„ is a two-component spinor of 
Grassmann variables which parameterizes the transfor- 
mation. This is a supersymmetry because the bosons and 



FIG. 5: The schematic phase diagram as a function of the 
ratio of the boson hopping tfj to the on-site boson repulsion 
energy U for a generic value of ho. 
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Although the evidences for the emergent supersymme- 
try, that is, the supersymmetric relation between cou- 
plings and the scaling dimensions, are obtained based on 
the calculation to leading order in e, we expect that the 
same conclusion holds to all orders in e, at least for small e 
for the following reason. The present model contains the 
supersymmetric Wess-Zumino theory in the sense that 
the supersymmetric Wess-Zumino theory can be reached 
at least by fine tunings of the bare parameters. There- 
fore, the Wess-Zumino theory should appear as a fixed 
point of this model in any case in the parameter region 
of Xi ^ ^ e although we don't know a priori whether 
it is a stable or unstable fixed point. Since the fixed 
point which was identified as the Wess-Zumino theory 
at leading order in e is the only fixed point which exists 
in that parameter range, and the higher-order terms in 
the e-expansion can not generate a new fixed point near 
the generic stable fixed point, the fixed point obtained 
to leading order in e should correspond to the Wess- 
Zumino fixed point to all orders in e. In other words, 
if the fixed point obtained to leading order in e were not 
the supersymmetric Wess-Zumino fixed point, another 
fixed point which corresponds to the Wess-Zumino the- 
ory should have appeared to leading order in e. 

Note that the supersymmetry transformations in Eq. 
(|12p mix bosons and fermions which carry different global 
U(l) charges, where the U(l) charges of the fermions 
and bosons are given by ~ 1 and — 2 respec- 
tively. Therefore the supercharges should carry the U(l) 
charge which is called R-charge. At the supersymmetric 
critical point, the super-Poincare symmetry is enlarged 
to an even bigger symmetry, that is, the superconformal 
symmetry which includes additional fermionic generators 
and the R-charge [1^. The additional fermionic genera- 
tors arise from the commutator of the supercharges and 
the conformal generators. The R-charge which enters in 
the superconformal algebra is related to the global U(l) 
charge as i? = — ^. The factor of 3 in the definition of 
the R-charge is due to the cubic superpotential of the 
Wess-Zumino theory (29|. Due to unitarity, there exists a 
constraint on the R-charge and scaling dimension of an 
operator. In 2+lD, the superconformal symmetry puts 
a lower bound on a scaling dimension as Do > \Ro\: 
where Dq and Rq are the scaling dimension and the 



R-charge of an operator O respectively. The equality is 
saturated for a chiral primary field and the one-loop cal- 
culation gives the exact anomalous dimensions for the 
fundamental fields with rjij, = r/^ — l/sfsO^I. Note that 
these coincide with the values obtained by putting e — 1 
in the one-loop results of Eq. (jl4p . Non-chiral primary 
fields generally receive radiative corrections and the crit- 
ical exponent for the order parameter is calculated to be 
j/^ = i + I + ©(e^) at leading order in ep[3|. 

IV. CONCLUSION 

In conclusion, we found a 2+ ID non-supersymmetric 
lattice model whose quantum critical point is described 
by the supersymmetric Wess-Zumino theory in an e- 
expansion. The supersymmetric critical point describes 
the generic second order phase transition between a nor- 
mal phase and a bose condensed phase of the bose-fermi 
mixed system. The exact anomalous scaling dimen- 
sions are predicted. In principle, the boson can dynami- 
cally arise as a Cooper pair in a system which has only 
fermions as microscopic degrees of freedom. In such a 
case, the critical point describes a superconducting phase 
transition. Therefore, supersymmetry can emerge at a 
critical point of a pure fermionic system. It is of interest 
to examine such possibility in the future. Finally, al- 
though the supersymmetric theory describes the generic 
second order phase transition in 4 — e dimension for a 
small e, in 2+lD wc can not exclude other possibilities. 
For example, the supersymmetric critical point may arise 
only as a multi-critical point due to an occurrence of 
other supersymmetry-breaking relevant perturbation, or 
the critical point itself may disappear due to a first order 
phase transition. 
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